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Abstract: We construct a basis set of infra-red and/or collinearly divergent scalar one- 
loop integrals and give analytic formulas, for tadpole, bubble, triangle and box integrals, 
regulating the divergences (ultra-violet, infra-red or collinear) by regularization m. D = 
4 — 2e dimensions. For scalar triangle integrals we give results for our basis set containing 
6 divergent integrals. For scalar box integrals we give results for our basis set containing 
16 divergent integrals. We provide analytic results for the 5 divergent box integrals in the 
basis set which are missing in the literature. Building on the work of van Oldenborgh, 
a general, publicly available code has been constructed, which calculates both finite and 
divergent one- loop integrals. The code returns the coefficients of l/e^,l/e^ and as 
complex numbers for an arbitrary tadpole, bubble, triangle or box integral. 
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1. Introduction 

The advent of the LHC requires a concerted effort to evaluate hard scattering processes 
at next-to-leading order in QCD. This requires both the calculation of tree graphs, for 
the leading order and real parton emission contributions, and the calculation of one-loop 
diagrams, for the virtual contributions. For most approaches to the calculation of one- 
loop amplitudes, the knowledge of scalar one-loop integrals is sufficient. For integrals 
with all massive internal lines these integrals are all known, both analytically ||l], |2|, |3| and 
numerically ^J. This paper therefore concentrates on integrals with some vanishing 
internal masses; these integrals can contain infra-red and collinear singularities. 

Many results have been presented in the literature for integrals with infra-red diver- 
gences regulated by introducing a small mass A for the divergent lines, especially in the 
important paper of Beenakker and Denner [^. Regulation with a small mass has been the 
method of choice in the calculation of electroweak processes in which massless lines are 
relatively rare. In QCD processes massless gluons are ubiquitous; in addition, light-quark 
lines are often treated as massless in the high energy limit. Consequently in QCD the 
method of choice for the regulation of collinear and infra-red divergent integrals is dimen- 
sional regularization. Therefore the results given in this paper are regulated dimensionally, 
although we shall exploit the relationship between the two methods of regulation where 
appropriate. 

Although some of the integral results presented in this paper are (to the best of our 
knowledge) new, the majority are not. However the analysis of the basis set required for 
a complete treatment of divergent box integrals at one loop is new. The results for the 
known integrals are dispersed through the literature and we believe it will be of use to 
collect the results in one place. Note that the general box integral in D-dimensions has 
been calculated in ref. [0]. However the results of that reference require considerable further 
manipulation before they can be used in practical calculations. 
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Given the complete set of tadpole, bubble, triangle and box integrals one can construct 
the scalar integrals for diagrams with greater than four legs |^, |l^. Thus a scalar 
pentagon in D dimensions, , can be written as a sum of the five box diagrams obtainable 
by removing one propagator if we neglect terms of order e 



i=l 



The general one loop A^-point integral m D = 4— 2e dimensions for > 6 can be recursively 
obtained as a linear combination of pentagon integrals |9| provided that the external 
momenta are restricted to four dimensions^. 



N 



i=l 



Thus for the purposes of next-to-leading order calculations, higher point functions A'^ > 4 
can be always reduced to sums of boxes. Note that for the case N > 7 the coefficients di 
in Eq. (|1.2| ) are not unique. 



2. Definitions and notation 
2.1 Definition of integrals 

We work in the Bjorken-Drell metric so that = ~ ~ ~ 'i- illustrated in Fig. || 
the definition of the integrals is as follows 



,4-D 



1 



nr'^rrJ (P — ml + ie) 

rD/2.2 2\ _ H- / rD, ^ 



I^{pf;mi,mi) = ^^ d^l 



in-rrJ (P - mf + ie){{l + qi)"^ - + ie) 

.A-D 

tD/ 2 2 2 2 2 2\ A* 

1 



X / d'^l 



{P — m\ + ie){{l + — 7712 + + (72)^ — m| + ie) 



tD / 2222 222 2\A'' 



X / d^l 



ITT 2 rr 
1 



{P -ml + ie){{l + qiY -m\+ ie){{l + - m| + ie){{l + q^^Y -ml + ie) ' 

(2.1) 

where g„ = Y17=iPi ^^^'^ 9o = and Sij = {pi +pjY- For the purposes of this paper we 
take the masses in the propagators to be real. Near four dimensions we use D = 4 — 2e. 



^Relations which have the same structure as Eq. (1.2) can also be derived, without the restriction that 



the external momenta lie in four dimensions. For the details of these relations we refer the reader to refs. 
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(For clarity the small imaginary part which fixes the analytic continuations is specified 
by +ie). /i is a scale introduced so that the integrals preserve their natural dimensions, 
despite excursions away from D = 4. We have removed the overall constant which occurs 
in D-dimensional integrals 



rr 



r2(l-6)r(l + e) l_ 

r(i-2e) r(i- 



1 - €7 + 



11 
2 



vr 
12 



O(e^). (2.2) 




Figure 1: The notation for the one- loop tadpole, bubble, triangle and box integrals. 

The final results are given in terms of logarithms and dilogarithms. The logarithm is 
defined to have a cut along the negative real axis. The rule for the logarithm of a product 
is 



ln(a6) = In a + In 6 + 77(0, b) , 

77(0,6) = 2TTi[0{-lm{a))0{-lm{b))9{lm{ab)) - e{lm{a))e{lm{b))0{-lm{ab))] 
So that 

ln(a6) = Ina + In 6, if Im(a) and Im(6) have different signs , 
In ^— ^ = Ina — In 6, if Im(a) and Im(6) have the same sign. 
The dilogarithm is defined as 

U2ix) = - / — In(l-z) 
Jo z 



1 s 
X X X 

12 + 22 + 



for Ixl < 1 . 



(2.3) 

(2.4) 

(2.5) 
(2.6) 
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A number of the most commonly useful dilogarithm identities are given in ref. |15]. In 
addition to the common identities, we also used an identity from ref. [p!^ ] in the calculation 
of the new box integrals 



Li2(l - ap2) + Li2(l - aQ^) - Li2(l - as) - Li2(l - at) 



U2[l 



+ Li2 1 



St 



+ Li2 1 



+ Li2 1 



(2.7) 



which holds for a = (P^ + q2 _ ^ _ t)/{P'^Q'^ - st). 
2.2 Analytic continuation 

The integrals given in sec. ^ are calculated in the spacelike region, Sij < 0,p? < 0. In this 
region, the denominator of the Feynman parameter integrals is positive definite and the i e 
prescription can be dropped. The analytic continuation is performed by restoring the ie 



Pi 



Pi +ie, 



Sij y Sij ~\~ 1 £ . 



rrii 



rrii — le . 



(2.8) 



The one-loop integrals which we present are all expressed in terms of two types of tran- 
scendental functions, 



Xi 



(2. 



ln|^JJa;ij , and Li2|^l-JJ; 

The individual Xj can be chosen such they that vary only on the first Riemann sheet of the 
logarithm, — tt < arg(xi) < tt. A complete specification of the continuation of terms of the 
form eq. (|2.9D has been given in ref. P]. The continuation prescription for the logarithm is 



ln( J|xi) ^Xl^"^' 



(2.10) 



\i=l 



i=l 



The continuation procedure for the dilogarithm is similar. For YVi=i -^d < 1 we use ||6|, |T2| 



Li2 1-n 



Xi 



i=l 



Li2 
6 



Xi 



i=l 



i=l 



In n^O -E^a 



\i=l 



i=l 



Li2 n^O i-ii^O 



(2.11) 



\i=l 



1=1 



1=1 



If 



nr=i 



> 1 it is expedient to make the transformation 



Li2 i-n^ 



-Lio 1 



i=l 



(2.12) 



We can then continue the resulting expression using eqs. (|1|, |2TI| ) as before. The con- 
tinuation procedure given in ref. [17] can be shown to be equivalent to the above. 
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3. Basis set of soft and collinear divergent integrals 



After Feynman parametrization and integration over we have for the triangle and box 
integrals 



tD/ 2 2 2 2 2 2\ 



/i^T(l + e) 



rr 



3 »i 

TT / dak 



1+e 



(3.1) 



tD/2222 -a a a a 



2 _2 



2 2^ /^'T(2 + e) 



4 .1 
TT / dafc 



where Y is the so-called modified Cayley matrix 

1 



rr 



Y 



+ m • - {qi-i - qj-iY 



2+e ' 

(3.2) 
(3.3) 



3.1 Landau conditions 



The necessary conditions for eqs. (3T,3^) to contain a singularity are due to Landau [18 
19]. If we introduce the bilinear form D derived from the modified Cayley matrix, 



D = '^aiajYij , 



(3.4) 



eqs. (|3.l| , |3.2| ) contain singularities if D = and one of the following conditions is satisfied 
for all values of j 

dD 



either a, =0 or 



dai 



0. 



(3.5) 



In general we have two classes of solutions of eq. ( |3.5D . In the first class we have 
solutions for the which are implicit functions of the masses and external momenta. Both 
physical and anomalous thresholds fall into this class. The leading Landau singularity 
occurs when = is satisfied for all j. Such solutions of the Landau equations will lead 
to divergences at the kinematic points where there are anomalous and physical thresholds. 

The second class of solution, which is of interest here, is the case where the external 
virtualities and internal masses have fixed values and the Landau conditions have solutions 
for arbitrary values of the other external invariants, Sij. Only these solutions will lead to 
soft and collinear divergences which are relevant for next-to- leading order calculations. 

As an example we consider the triangle shown in Fig. (|2|a) which contains a soft 
singularity. In this case the denominator is given by 

D = {ml + - pl)a2a3 + mlal + m|a| . (3.6) 

This expression satisfies the Landau conditions for 02 = as = and ai arbitrary. A second 
example is the triangle shown in Fig. (||b) which contains a collinear singularity. In this 
case the denominator reads 



D = {ml- ^2)0203 + {ml - pl)aias + m|a| , 
which satisfies the Landau conditions for 03 = and ai, 02 arbitrary. 



(3.7) 
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3.2 Soft and collinear divergences 

Prom the Landau conditions it follows that a necessary condition for a soft or collinear 
singularity is that for at least one value of the index i [pH] 



Yi+i ^+l = Yi+i i+2 = Yi+i 1 = 0, soft singularity , 
Yi i = Yi+i j+i = Yii+i=0 , collinear singularity . 



(3.8) 
(3.9) 

The indices in eqs. ( p.8| , p. 9] ) should be interpreted mod A^, where N is the number of 
external legs. Thus the structure of the Cayley matrices for integrals having a soft or 
collinear divergence is as follows 



soft 



/... ... 

.. ... 



\ 



Y. 



collinear 



7 




(3.10) 



In order to have a divergence, we must have at least one internal mass equal to zero, i.e. 
at least one vanishing diagonal element of Y. 

3.3 Basis set for triangle integrals 

We list here the basis set of six divergent triangle integrals as shown in Fig. |3|. All other 
divergent triangle integrals can be derived from this set. First we have two integrals with 
no internal masses, 

1. /3^(0,0,pi; 0,0,0) 

2. /3^(0,pi,pi,0; 0,0,0). 

Second, we have three integrals with one massive internal line, 

3. /3^(0,p2,p2.o^o,m2) 

4. /3^(0,p2,r7i2;0,0,m2) 

5. If (0,m2,m2;0,0,m2). 
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Figure 3: The six divergent triangle integrals. Lines with a zero internal mass mi = (for internal 
lines) or a zero virtuality, = 0, (for external lines) are shown dashed and are unlabelled. Solid 
lines have a non-zero internal mass, or a non-zero virtuality. Lines with the same color have the 
same internal mass and/or virtuality. 



Last, we have one integral with two massive internal lines, 

6. If (mi, pi, mi; 0, mi, mi). 

The set is in fact overcomplete since the second integral can be obtained from the m^ — ?■ 
limit of the third integral. However from a numerical point of view it is expedient to 
categorize the integrals by the number of vanishing internal masses and to treat the two 
cases separately. 



3.4 Basis set for box integrals 

In this section we demonstrate that a basis set sufficient to describe all box integrals 
with collinear or soft singularities can be constructed from sixteen integrals, illustrated in 
Fig. All other divergent box integrals can be derived from this set. ^ The integrals are 
characterized by the number of internal masses which are equal to zero. Each one of these 
divergent integrals has a characteristic modified Cayley determinant as shown in Fig. ^ 



satisfying the conditions of eqs. (3^, 3.9 ) 



As discussed in sec. 3.1, at specific kinematic points there are threshold singularities derivable from the 



Landau conditions which can lead to singular integrals not derivable from our basis set. 
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Figure 4: The sixteen divergent box integrals. Lines with a zero internal mass nii — (for internal 
lines) or a zero virtuality, p| = 0, (for external lines) are shown dashed and are unlabelled. Solid 
lines have a non-zero internal mass, or a non-zero virtuality. Lines with the same color have the 
same internal mass or virtuality. 

3.4.1 Integrals with no internal masses 

There are five integrals with no internal masses, 

L If (0,0, 0,0; si2,S23; 0,0, 0,0) 

2. /f (0,0,0,p|;si2,S23;0,0,0,0) 

3. If (0,p2,o,p2 ;s^2,S23; 0,0, 0,0) 

4. If (0,0,pi,pl;si2,S23;0,0,0,0) 

5. I^iO,plplpl;si2,S23;0,0,0,0). 
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Figure 5: The structure of the modified Cayley determinant for the sixteen divergent box integrals. 

This set of integrals is sufficient to describe all the integrals with massless internal lines, 
but it is overcomplete since in ref. an expression is given for box integral 3, which can 
also give results for integrals 1 and 2 by taking the appropriate limit. However, as before, 
for reasons of numerical expediency it is convenient to retain the overcomplete basis. 



3.4.2 Integrals with one non-zero internal mass 

If we have one non-zero internal mass we can take this without loss of generality to be the 



last one (7714). In this case the modified Cayley matrix, eq. (3^) is 
/ 



Y 





2 1^2 1 



2P1 


2 1 , 







4 - 2P4 2"^4 - 2*23 2"^4 



i„2 
2P3 



1 2 

1 2 
1 2 



1*23 
2P3 



(3.11) 



With si2,S23 fixed we can apply four conditions to potentially create a soft or collinear 
divergence, namely 

Pi = 0, P2 = 0> pI = w-i pI = 'ml- (3.12) 

However performing the interchange p\ -H- ^ p1, (with fixed) corresponds to a 

relabelling of the diagram. In addition setting either |?| = m| without setting p| = 0, or 
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pI = ml without setting Pi = does not lead to a divergence. If we denote the apphcation 
of the four conditions, eq. ( 3.12| ) on Pi,P2,P3,p1 by {i,j,k,l) we have following 15 cases: 



6. ) 

7. ) 



10.) 



(1,2,3,4) 

(1,2,3)^(1,2,4) 

(1,2) 

(1,4) = (2, 3) = (1,3, 4) = (2, 3, 4) 
(1) = (2) = (1,3) = (2,4) 
(3) = (4) = (3, 4) = finite . 



(3.13) 



To be quite explicit, the notation (1, 2, 3) corresponds to setting = 0,P2 = andp| 
etc. The first five lines in eq. ( |3.13| ) correspond to the integrals 6-10 given below, 

6. 4^=^"^'^(0, 0, m2, m^; su, S23; 0, 0, 0, m^) 

7. /i^='-''^(0, 0, m\pl; su, S23; 0, 0, 0, m^) 

8. 4^='-'^^(0, 0,pIpI; si2, S23; 0, 0, 0, m^) 

9. iI''=^~^'\0,pIpI m2; su, S23; 0, 0, 0, m^) 
10. iI''=^~^'\0,pIpIpI; si2, S23; 0, 0, 0, m^). 

The last case in eq. ( p.l3D does not lead to a divergent integral. 



m 



'4! 



3.4.3 Integrals vi^ith two adjacent internal masses 

Without loss of generality we can take the two non-zero adjacent internal masses to be 7713 
and 1714. In this case the modified Cayley matrix is 



Y 











1^2 1^2 



12 12 

2 777-4 — 2P4 



\ 



^mi-isi2 ^mi-ipi 7ni ^Tni + i7n| - 

V Irnj - IpI \ni\ - \s23 \ml + \vn\ - \pl m\ 



(3.14) 



A necessary condition to have any divergence is p\ = 0. This gives the integral 13. Applying 
either p\ = m\ or p\ = m\ gives a pair of integrals related by relabelling, integral 12. 
Applying both p^ = m\ and p| = 777^ gives integral 11, 



11. 4^ ^ ^'^(0,r7l|,p|,77l|;Si2,S23;0,0 



,7773,777,4; 



12. /, 



{D=A-2e}, 



(0, ?77|,P3,P4; S12, S23; 0, 0, mi, mi) 



13. 4"""^ ^'\0,pIpIpIsi2,S23;0,0 



,m^,ml) 
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3.4.4 Integrals with two opposite internal masses 

Without loss of generality we can take the two non-zero opposite internal masses to be 7712 
and m4. In this case the modified Cay ley matrix is 



Y 



/ n 1^2 i„2 i„ 1^2 i„2 \ 

^mj - \pI ml \ml - \pl \ml + \ml - \s2s 



(3.15) 



Here we can only have a soft divergence since there is no pair of adjacent zero internal 
masses. Setting p\ = ml, pi = m\ or p^ = rnl^p^ = m\ gives two integrals related by 
relabelling (15). Setting both conditions gives integral 14, 

14. /j'^~^~^^-^(m^, m^, ml, m|; S12, S23; 0, m^, 0, mf) 

15. li'^'^''^^^ {rnl,pl,pl,ml-, S12, S2'i,;^.,ml,{),ml) . 
3.4.5 Integral with three internal masses 

Without loss of generality we can take mi = 0. With only one zero mass, there can only 

be a soft singularity. This requires the two adjacent external lines to satisfy the conditions, 

2 22 2 
p{ = m2,pi = m| 

16. /f (m|,p|,p|,m|;si2,S23;0,m|,m|,m|) . 

The concludes the listing of our basis set of the divergent box integrals. 

4. Results for integrals 

4.1 Tadpole integrals 

The result for the tadpole integral is given by 

/f (m^) = -^2e p(_^ ^ j^2 _ .^jl-e 

The 1/e pole corresponds to an ultraviolet divergence and the analytic continuation has 
been made explicit in this case. 

4.2 Bubble integrals 

Following 't Hooft and Veltman [|l|, the result for this bubble integral in our notation is, 
rD/-„.2 2\ _ /^^''r(e) 2 I /I „\_2 



h (s; rni,m2) = / [-7(1 - 7)5 + 7m2 + (1 - 7)"ii - i^] ' 

rr Jo 

= ^^^'{^~ ln(-7(l-7)s + 7m^ + (l-7)m^-ie)| +0{e) 



- 12 - 



.2e 



e Jo ^ s s 



.2e 



+ 2- ln(s - ie) + ^ [ 



7.1n(^)-ln(7.-l] 



i=l 



where 71 2 are the two roots of the quadratic equation, 



m: 



71,2 



+ m1± y^(s — ml + mfy — 4s(m^ 



2s 



The special hmits for this integral are 

r 1 



/f(s;0,m^) 



In 



If (s; 0,0) 



+ 2 + 



' 'i + 2!>+0(e) 



+ 0(e) 



+ 0(e) 



(4.2) 



(4.3) 



(4.4) 



As before, the 1/e pole corresponds to an ultraviolet divergence and the analytic continu- 
ation has been made explicit in these cases. 



4.3 Divergent triangle integrals 

In this section we give the explicit results for the six divergent triangles. The results 
for the triangles have been presented already by many authors and are given here only 
for completeness. The results are reported in the spacelike region below all thresholds. 
The analytic continuation is performed using the prescription given in sec. 2.2. Each 
expression for a triangle integral stands for the 6 different labellings of the triangle obtained 
by repeated application of the following identities 

tD/ 222 2 2 2\ rD/ 222 2 2 2\ 

h (Pl>P2>P3;"il'"i2'"i3) = ^3 (P2'P3'Pi;"^2'"i3'"^l) ' 
tD/ 2 2 2 2 2 2\ tD/ 2 2 2 2 2 2\ 



4.3.1 Triangle 1: {0,0, pj; 0,0,0) 



/3^(0,0,p^; 0,0,0) = ^ 



pI 



-pI 



+ -In^ 



-pI 



4.3.2 Triangle 2: ll^ {0,pi,pi;0,0,0) 



+ Oie) 



(4.5) 



2e 



I^{0,pIpI,0;0,0,0) = ^ 



P2-pI\^ 



-A 

-pI 



1 

+ 2 



In^ 



pI)-' 


-i-pl 


^1 


pI 


-pI 










-pI\ 


ln2(- 


-pI\\ 


ii"^ ) 


^2 )\ 



+ 0{e) 



(4.6) 
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In the limit ~^ P3 define r = (p| — pI)/?*! ^'^'^ following expansion 



^3 (0,P2,P3,0; 0,0,0) - 2 



'^l-;Ulnf^U^fl + ln^-^^ 



+ 0(e,r2). 
(4.7) 



4.3.3 Triangle 3: /f (0,p|,p|; 0, 0, m^) 



2 \ f 



iln(4^ 
e V — P2 



+ Li.(4) 



+ In 



2 2 
2 ("rn/^-j4 

9 



In^ 



2 2 
m — P3 



(4.8) 



In the limit ~^ P3 define r = (p| — P2)/("^^ ~ P2) ^^"^ following expansion 
/3^(0,pi,p^,0; 0,0,0) 



m? — P2 



1 - - - - In ^ 



p2 



In 



r 
2p: 



- 2plm^ - pi (m^ -pI\ , 2 , 2 
^ In I I +m +P2 



(4.9) 



Rewriting eq. ( f4.8D in the following form makes the m — )■ limit and the agreement 
with eq. (W^ manifest 



l^{0,plpl;0,0,m') 



1 



P2 — P3 I e ^ JTi^ — P2 ^ 2 - 



+ In 



(m^ - p2) ^^ ^ -p2(^2 _ pi) ^ ^ (^i _ pi) 



-pi 



In 



In 



,2 ; V „2 

pi{m^ -pI)' 



m' 



Li2(— ]+Li2(^)KO(6). 



(4.10) 



4.3.4 Triangle 4: /f (0,p|, m^; 0, 0, 



m 



/3^(0,pi,m2;0,0,m^) = ( il^ / ^ 2 



2 \ ^ 



1 



1 1 , 

^ + ~ I 2 2 
ze^ e V — P2 



^ 1,2 

H h - In^ 

12 2 



rrfi — P2 



Li2 



-pi 



?7l2 — P2 



4.3.5 Triangle 5: if (0, m2, m2; 0, 0, m2) 



+ 0{e) 



(4.11) 



73^(0, m2,m2;0,0,m2) 



1 



rn-2 / m2 \ 2e 



+ 1 



(4.12) 
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4.3.6 Triangle 6: {ml, s,ml;0,ml,ml) 

The result for this triangle integral can be obtained from ref. eq. (C3), by the normal 
replacement rule |22] which is true in the case of a soft singularity 

InA^ ^ ^ + ln/i2 + 0(e). (4.I3) 

For s — (m2 — m^)'^ 7^ we have 

q '{m2,s,m^;0,m2,m2) = 

m2m3(l — Xg) 

I r 1 1 5 

x<ln{xs) -ln(j;s) +21n(l -x^) +ln(- 




- ^ + Li2(x2) + i ln2 + Li2(l - xs^) + Li2(l - x,^) I + 0{e) , (4.14) 
2 m3 1712 I 

where Xg = —K[s + ie, m2,rn^) and is given by 

K(z,m,m) = — \ ^ z^{m-mY 

K{z,m,m') =—1 z = (m — m')^ . (4-15) 

For s — (m2 — m-a)^ = this becomes 

{D=4}/ 2 2 n 2 2\ 1 

3 V 2> ' 3' ' 2> 3; 2m2m3 

f - 2 - !^^i±^ In f ^1 I + . (4.16) 

4.4 Divergent box integrals 

In this section we give the explicit results for the sixteen divergent boxes. The results 
are reported in the spacelike region below all thresholds. The analytic continuation is 
performed using the prescription given in sec. 2.2. Each expression for a box integral 
stands for the 8 different labellings of the box obtained by repeated application of the 
following identities 

tD/ 2222 222 2n tD/ 2222 222 2\ 

h (Pl,P2>P3'P4.Sl2,S23,mi,"i2,"l3'"^4) = H (P2 ) ^3' P4> Pi > •523 , •§12 , ^2, mg , 711,4, ) > 
tD/ 2222 222 2\ tD / 2222 222 2\ 

h (Pi,P2)P3'P4) •512, •S23,"^l, 7712, mg, 711,4) = h [Pa, Ps^ P2^ Pv ^u, S23, m^, ui^^, m^, 1712) . 

(4.17) 

Where we have found the integrals in the literature we give references. To the best of our 
knowledge the results for boxes 9, 10, 11, 12, 13 are new. 

4.4.1 Box 1: 4^(0, 0,0,0;si2,S23; 0,0, 0,0) 

„2e 

/4^(0,0,0,0;si2,S23;0,0,0,0) - 



•S12.S23 



-S12)-' + (-523-)-^ ) - In' (^l^) - K 0{e) . (4.18) 



This result is taken from 111 
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4.4.2 Box 2: if (0, 0, 0,p^; S12, S23; 0, 0, 0, 0) 



2. 



If (0,0,0, p^;si2,S23; 0,0, 0,0) 



2e 



S12S23 



||((-512)-^ + (-523)-^ - i-pl)-') - 2Li2(l - ^) - 2Li2(l 



S23' 



In^ 



Sl2\ TT 



-S23 



(4.19) 



This integral is given in p3| , ^]. An alternative formulation with three dilogarithms is 
given in 



4.4.3 Box 3: If (0,p^, 0,pl; S12, S23; 0, 0, 0, 0) 



if (0,pi,0,pl;si2,S23;0,0,0,0) 



2€ 



S23'S12 - pIpI 



{I ((-^12)-^ + (-523)-^ - {-pD- - i-pl)-') 



2U2 1 



+ 2Li2 1 



i4_ 
pIpI 

S12S23 



2Li2 ( 1 
-In 



P2 



S23 

2 I -■512 



2Li2 1 



S12 



2Li2 1 



P4_ 

S23 



-S23 



+ 0(e). 



(4.20) 



This result is taken from [11]. As for all the integrals, the analytic continuation of this 



result follows the procedure detailed in sec. 2.2. An alternative form in which the analytic 
continuation is manifest is given in ref. pi]]. 



When the denominator in the overall factor in eq. ( 4.20 ) becomes small and the entries 



in the two pairs {su, S23) and {P2,pi) have the opposite sign we can expand in r = 1- 



, P2P4 

S12S23 



where 



If (0,p^,0,p^;si2,S23; 0,0, 0,0) 



S23S12 



\2 + r) + (2-^r) + (2 + r)(ln, „ 
e 2 \ \ fi^ 



-Pa 



+ 2 



S23 



Sl2 



+ r 



S23 



Sl2 



Uiz) 



ln(z) 



U(z) 



Ln (z) + 1 



+ 0(e,r2), (4.21) 



(4.22) 



1-z' 1-z 

Thus in this region the residue of the overall pole at S23S12 = P2P'i vanishes and we obtain 
a numerically stable expression. The other region si2,S23 > and p^'Pl < '^i^^ 
versa (512,^23 < and > 0) is the region of the Landau pole which gives a large 

contribution to the imaginary part. 
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4.4.4 Box 4: It' {0,0, p^, pi; su, S23; 0,0, 0,0) 

,2e 



I^{0,0,pi,pi;si2,S23;0,0,0,0) 



„2 2 



■((-^12)-^ + (-^23)-^ - {-pD- - {-pId + ^{{-pD-\-plr')i{-si2r 



2Li2 1 



pI 



2Li2(l-^)-ln2 

S23/ V-'523 



-Sl2 



pi_ 

S23 



This result is taken from ref. [TJ|. (See also ref. pi|). 
4.4.5 Box 5: I^{0,pl,pl,pl,si2,S2z;0,0,0,0) 

,2e 



+ 0{e). 



(4.23) 



I?{0,pi,pi,pi;si2,s2z;0,0,0,0) 



(S23S12 - pIpI) 

X ||((-si2)-^ + {-S2,r - {-pir - {-pir - i-pi)-') 
^((-P2)~^(-pi)'^)/(-«23)-^ + ^^{{-plr{-plr)/{-si2)- 



2U2 1 1 



(4.24) 



P2_ 

S12 



2U2 1 



P4_ 

S23 



+ 2Li2 1 



phi 
S12S23 



In^ 



-■S12 

-S23 



+ 0(e) 



This result is taken from [11]. As for all the integrals, the analytic continuation of this 
result follows the procedure detailed in sec. 2^. An alternative form in which the analytic 
continuation is manifest is given in ref. [pl| ]. 

When the denominator in the overall factor in eq. (4.24) becomes small and the entries 

2 2 

in the two pairs (si2, S23) and (pIjPI) have the opposite sign we can expand in r = 1 — j^2£4_ 



lFiO,pipi,pt si2,S23;0,0,0,0) 



,i2 ^2 J2. 



In 



S23S12 
2 



-S12 



In 



-Ps 
-S23 



2-(l + ^)Lof^ 

S23 V'S23 



+ r 



P4_ 

S23 



Ln 



A. 

S23 



+ 0(e,r2). 



(4.25) 



with Lo,Li as in eq. ( 4.22 ). 

4.4.6 Box 6: if (0, 0, m^, m^; S12, S23; 0, 0, 0, m^) 



{D=4-2e}/„ n 2 2 n n n 2 

i4 '(0, 0,m ,m ;si2,S23;0, 0,0, m 



1 



2 \ e 



e2 e 



l(21n( 



) + ln(- 



Si2(m.2 - S23) V"^ 
+ 21n( 2 >H 2") - 



+ Oie) . 
(4.26) 



The result for the real part of this integral in the region S12 > 0, S23 < is given in ref. 
eq. (A4) (note differing definition of e). 
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4.4.7 Box 7: /^^(O, 0, m^p|; si2, S23; 0, 0, 0, m 

,2 



If (0,0, m^p^;Sl2,S23; 0,0,0, m") = ( ^ 



,2 2 



3 J_ _ 1 

2? ~ 7 



21n( 1-^ ) + ln 



m^y S12(S23 — "i^) 



«|)-ln=('l-4 



5^ 
T2" 



(4.27) 



This integral is obtained from eq. ( A4) (first equation) of ref . |25| ] . (The real part was given 
earlier in eq. (6.75) of [p^].) 



4.4.8 Box 8: if (0, 0,p^,p^; S12, S23; 0, 0, 0, m 



.2 Ji. 



If (0,0,p^,p^;si2, S23; 0,0,0, m^) 



.2 2 



1 



1 1 



-2 Li2 1 



In — ^ + In 



S12(S23 - m^) 
(m^ - S23)^ 



2 2 
m — P3 

- S23 



^^lj^2/--12 



2 Li2 1 



6 2 
-In 



-In^ 



-S12 



2 2 
m — p| 

- S23 
+ 2 In 



Si2?7l2 



-•S12 



2 2 

m — P3 



In 



2 2 
m — P3 



ln(!^)ln 



In 



- S23 



2 2 
m — 



+ 0(6). 



(4.28) 



This integral was constructed from the expression eq. (B6) of ref. [p7|], which is valid 
for the real part in the region S12 > 0, S23 < 0. 



4.4.9 Box 9: /f (0,pi,p|, m^; S12, S23; 0, 0, 0, 



+ Li; 



(0,pi,pi m^; si2, S23; 0, 0, 0, m- 
^ — ^3)(m^ 



1 1 f S12 {m? - S23) 
In 



fi2(s23 - "T-^) 1 2e2 e Vp2 
V P2 ^ 12 V /^"T- ^ J 



(4.29) 



4.4.10 Box 10: if (0,p2,p2^p2. g^^, S23; 0, 0, 0, m2) 



h (0,P2)P3'P4;'5l2,S23;0,0,0,m 
,2 „2\„2 



X 



(S12S23 - m'^su - pIpI + rrflpl) 
r 1 / + Li2 (1 + -^^3)A _ Li + (rn^^plKm^x 

[e V (?n^ — S23)si2^ ^ ^2"^"^ / V Si2m^ / 



+ 2Li2(l - ^^^) - 2Li2(l - i^) + 2Li2(l - 

,2 ^2\^2 

/ I I.I 1 1. \ /{II 

+ 2 In 



)'"(T!S^)h°(')' 



. - S23 ^ ^ ("^^ - S23)S12 



(4.30) 
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4.4.11 Box 11: I4 {0, ml,pl, ml; si2,S23;0,0, ml, ml 
If (0,mi,pi, ml; .12,. 23; 0,0, mi, ml) 



(m| - si2)(m| - S23) 



1 1 /(m|-S23)( 



mt 



S12) 



m^m^jw^ 



vr 



+ In 



2 

Km^J 



1 



73I 



734 



2 In 



■In^ 



mn 



S12 



m^fi 

73^4 



In 



ml - S23 



m^fj, 



734 - 1 



Oie) , 



where 



2 2 
m^ — mj 

— ± 

2 



'(1 



m^-^m^ 
P3 



4m| 



and + 7^^ = 1. Assuming ml > mg, in the hmit P3 — )• we obtain. 



73+4 



ml 



734-1 ■ " 734-1 "^3 

and this expression reduces to the form given in eq. (6.77) of Hopker 

1 



/4 (0,m|,0,m|;si2,S23;0,0,mi,ml) = — ^ 



(mi-.i2)K-S23) 



1 1 



■In 



+ 21n 



m 



K-S23)H 



S12 



S12 



m^fi 



In 



2 

'4 — ^^23 

m^fj, 



m^m^fi 

2 

m 



4.4.12 Box 12: /f (0, mi,pi,pl; S12, S23; 0, 0, m^, m^ 



(4.31) 



(4.32) 



(4.33) 



(4.34) 



if (0, m|,p|,pl; S12, S23; 0, 0, m|, ml) 



(^12 -mi) (.23 -m|) 



1 
2? 



(^ (ml-.23)(rni-.i2) ^ , .,^(mt-S2s 

msfi 



In 



2 / ml - p1 



ms^u 



{ml - pl)m,ii 
± (ml-vl 
12 



In 



mi - S12 



ma/i 



m% 



S\2 



In 



m. 



m; 



7^4 



2 V73+-I 



■In^ 



734 



734 
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2Li2 1 



ml 



(ml - S23)^ 



Lio 1 



(ml - p1) 7+ 
(mi - S12) 7I3 - 1 



Lio 1 



where 7^^ is given in eq. ( g32|) . In the hmit pi = 0, we get 

/f (0,mi,0, p1; si2,S23; 0,0, mi, ml) = 7^ 2T 

(Si2 — m3j(S23 — m4j 



(ml - pD 743 \ \ 
(mi - S12) 7^3-1// 
(4.35) 



2^-7^" 



(ml-S23)(mi-.i2)^^21n^"'^-^23 



(ml-pl)m3M 



m3/i 



In 



m3 - S12 
m-sfi 



-19- 



.„M^Vg..n(^)l.(;|)-l.„^(| 



- 2 Li,(l - i^^^) - Li,(l - f!^^) - Li,(l - ^p^fLEiL) I + 0(.) , 

(4.36) 

4.4.13 Box 13: if (0,pi, pi, pi; S12, S23; 0, 0, m|, m|) 

/f(0,pi,pi,pl;.i2,.23;0,0,mi,ml) = l|iln H " l^i)K - p1) 

A[e V("i3 - ■Sl2)("^4 - «23) 

H (mi-si2)/ (?"I-S23) 73^4-1^ ("^I-S23)7:r4-1^ 

- 2UJ1 - j4^) - LiJl - ) - Li2(l - ^f4^^^^) 

V (m|-S23)^ ^ ("^3 -512)723 - 1^ ^ ("^3 -^12)743 - 1/ 

+ 2Li2(i - ^ 2in m 

V (m§ - si2)(m| - S23)^ V /i^ / V /i^ / 

_ 1^2 1^1 \ _ 1^2 ^ f rnl-pl ^ /mfx ^ / m^-p^N /m|x 

\ jl^ J V /X^ / ^?TT,| — S23^ V ^2 y Vm| — Si2^ V ^2 / 

- ^ In^ (^) - ^ In^ (^^) I + , (4.37) 
2 ^734 -1^ 2 V^3^-l/J 

where 7^^ is given in eq. ( [4.32 ) and 



A = (S12S23 - "1-^523 - "T'lsi2 - vIpI + '^IpI + "^4^2) 

= {ml - si2){ml - S23) - {ml - pI){'^1 - pI) ■ (4-38) 



In the hmit P3 — >• this simphfies to 

1 fli ( {ml-pl){ml-pl) 



-^f (0, pi, 0,p|;si2,S23; 0,0, m|,m|) = ^ j^ln 



(mi-.i2)(ml-.23), 

2 Li2 (1 - ^4^) - (1 - - (1 - i 

V (m^-si2)/ V (m|-S23)^ ^ rn|(m|-S23)^ 



2Lio 1 ^ 



2 



^) - Li2 (1 - _ Li2 (1 - 4-^4^) 



{ml - S23) ^ ^ {ml - S12) / V m| (m§ - S12) / 

+ 2Li2(l - ) + 21n In 

V (m§ - si2)(m| - S23)^ V /i^ / V 

- In^ f - In^ f + In (4^) 1^ (4) + 1- (4^) In (4) 
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4.4.14 Box 14: /f (m^, m^, m|, m|; S12, S23; 0, m|, 0, m|) 

We can obtain this doubly IR divergent box integral from eq. (2.13) of ref. H, using the 



simple replacement rule in eq. (4.13). We obtain 

rD. 2 _2 _2 „ .n_2n_2N_ "2 2:23 ln(x23) [ 1 



tU I I a a I nZrv2\ 

^4 [m2,m2,m^,m^; S12, S23; 0, mg, 0, m^) 



m2miSi2 1 — 



23 



+ In 
e V-S12 



The variable X23 is defined in terms of the function eq. (4.15), such that 



X23 = --f^(s23 + m2,m4) . 
In the limit S23 — (^2 — rn^)'^ — )• we have 



(4.40) 
(4.41) 



{ml, ml, ml, ml; S12, S23; 0, ml,0, m; 



m2miSi2 I e 



+ In 



-S12 



+ 0(e,(l-X23)'). 

(4.42) 



4.4.15 Box 15: /f*(m^,p^,p|, m|; S12, S23; 0, m|, 



We can obtain this IR divergent box integral from eq. (2.11) of ref. Q, using the simple 
replacement rule eq. ( [4.13 ). We obtain 



h ("i2>P2;P3) "^4; •S12, S23; 0, 0, 7714) 

1 1 



2:23 



7n2m4Si2(l — x"^ 



X < lnX23 



Li2(l - X 



e 2 



In 



2^23 - In (— — ) - In ( 



237 

2 2 2 2 



-S12 



-S12 



23, 



+ ^ln''y+J2 Li2(l - ^23^0 I + 0{e) 
^ P=±i ) 



where 



7772 (ml - pI) 
7774 (7^2 — P2) 



(4.43) 



(4.44) 



and the variable X23 is defined in eq. ([4.411) . 

icf 

T-U/ 2222 n2rv2 
h ("^2'?'2)P3' S12, S23; 0, 7772, 0, ^4 



For 7774 — P3 small it is useful for numerical purposes to rewrite eq. (4.43) in the form 

J2 -2 ^2. . „ . r, „„2 r, „„2\ ^23 



"727774512(1 - ^33) 



X < lnX23 



i-lnx23-ln(4)-21n(!^i^ 



Li2(l - x^a) + Li2(l - X23y) - Li2(l -—))+ 0{e), 

X2?. 



(4.45) 



and similarly for 7773 — pi small. 
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In the limit X23 — >• 1 (i.e S23 = (m2 — ^4)^) we obtain 

1 



tDi 2222 n n 

h {'m2,P2,P3^ "^4; si2, S23; 0, m2, 0, 771,4) 



27772 7774512 



X I i + In (^) + In r'^-Ph + In f - 2 - -i±^ InJ 

I e V77727774/ V — si2 / V — S12 ^ (1 ~ y) J 

+0(e,(l-X23)'). (4.46) 
4.4.16 Box 16: /f (r77|,p|,p|, 777^; S12, S23; 0, 777^, 777^, 777^) 

We can calculate this IR divergent box integral from eq. (2.9) of ref. using the simple 



replacement rule eq. ( [4.13 ). We obtain 



tD/ 9 9 9 9 9 9 9\ "^23 

h (m'2;P2'P3'"i4;Sl2,S23; 0,7772, 7773, 7774) = ^ 

77727774(512 - "T'3)(l " ^23) 

X ( - - 21n(x23) In (^^) + In2(x2) + In2(x3) - Li2(l - xl,) 

I e \r77|-si2/ 

, T - n \ , T - f-i ^23 X , T - n a;23a;2x , j ■ a^232;3J , ^ 
+ Li2(l - X23X2X3) + Li2(l ) + Li2(l ) + Li2(l ) > + 0{e) 

X2X3 X3 X2 



(4.47) 



where X23 = --f^(s23, "1-2, "1-4), X2 = -K{p\,m2,m^) and X3 = -i^(p|, 7)13, 7)1,4). 
In the limit X23 — )• 1 (i.e S23 = (7772 — 7774)^) we obtain 

tD/ 2222 n22 2\ 1 
h ("72;P2;P3; "^4; •512,523; 0,771,2, 7773, 7774) = ^ 



1 , / 7773yU \ 1 + X2X3 

X <^ - + 2 In ' ^ ^ ^ 



2r7727774(si2 — 7773) 

^ ^2 



ln(x2) + ln(x3) - ln(x2) - ln(x3) 

X3 - X2 



.r77|-Sl2^ 1 - X2X3 

+ 0(e,(l-X23)'). (4.48) 

Special choices of p^iPi ^nd (non-zero) values of the masses m^^m^^rn^ will not lead to 
further divergences. 

4.5 Special cases for box integrals 

In this section we give some examples of non-singular limits of the basis set of box integrals. 
The first example illustrates that the basis is overcomplete. If we look at the 777^ — )• limit 
of box 8, eq. ( 4.28| ), we find after little work that it reproduces the result for box integral 4, 



eq. ( [4.23| ). That this should be the case is clear from the form of the modified Cayley matrix 
y, as shown in fig. ^ Taking the limit m? — )• does not introduce any new singularities 
of the form given in eqs. (3.8, |3.9|). In a similar way one can show, for instance, that box 



13 reduces to box 10 in the limit 7773 — t- and box 10 goes to box 5 in the limit 7774 — )■ 0. 
Since these limits, which are analytically simple, can sometimes be numerically delicate, 
we choose to treat these integrals as separate cases and to categorize the integrals by the 
number of the internal masses which are non-vanishing. 
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Divergent box 



Special Case 



Reference 



Box 8 



^4 



0,0,pi,p|;si2,S23;0,0,0,m2) 
0,0,0,p|;si2,S23;0,0,0,m2) 



m, eq. (6.71) 



m, eq. (A17) 



Box 11 



^4 
-'4 
-'4 



0, m|, 0, m|; si2, S23; 0, 0, m|, m|) 
0, m^, 0, m^; S12, S23; 0, 0, m^, m^) 
0, ni^^pl.m?; S12, S23, 0, 0, m^, m^) 



126], eq. (6.77) 



l2i, eq. (6.70) 
[H, eq. (30) 



|2§], eq. (A4), third eqn. 



0, eq. (B7) 
m, eq. (6.74) 
eq. (6.78) 
25], eq. (A4), second eqn. 
25], eq. (A4), fourth eqn. 



Box 12 



-'4 
-'4 
-'4 

-'4 



0, m|, 0,p|; S12, S23; 0, 0, ml, ml) 
0, m|, 0, m|; S12, S23; 0, 0, ml,ml) 
0, m|, 0,p|; S12, S23; 0, 0, m|, m|) 
0, m^, 0,pI; si2, S23; 0, 0, m^,m^) 



0,m^,p^,p|;si2,S23,0,0,m^ 



Box 13 



-'4 
-'4 

^4 



0,p^ 0,p2; S23; 0, 0, m\m^ 



0, m|, 0, m§; S12, S23; 0, 0, m|, m|) 
0,0,p|,pl;si2,S23;0,0,?n2,m2) 



^, eq. (6.72) 
I eq. (6.79) 
1 (v2), eq. (A19) 



Box 16 



-'4 
-'4 
-'4 
-'4 
-'4 



m'^,pl,pl,m'^]Si2,S23]0,rn'^,rn'^,'rn^) 
m^, 0,p|, m^; S12, S23; 0, m^,m^,m^) 
m^, 0, 0, m^; S12, S23; 0, m'^,m'^, iv?) 
ml, 0, 0, mj; S12, S23; 0, ml, ml, ml) 
ml, 0, 0, S12, S23; 0, ml, ml, mj) 



25], eq. (A4), sixth eqn. 
H, eq. (A4), fifth eqn. 
H, eq. (A3) 



\2§, eq. (6.73) 



\2§, eq. (6.76) 



Table 1: Special cases of the 16 basis integrals available in the literature. 



The second illustration gives an example of an integral which is obtainable from one of 
our basis set integrals by taking a non-singular limit. If we look at the first entry of tab. || 
and take the limit p^ = p| we reproduce the result of ref. 126], eq. (6.71). A little care is 
required since ref. ]26] only gives the result for the real part of the integral in a physical 
region and has a different e-dependent overall factor. Table || details the examples which 
we have found in the literature which are non-singular limits of our basis integrals. We 
have checked that we are in agreement with all these special cases. 



5. Numerical procedure and checks 

We have constructed a numerical code which for any A^-point integral returns the three 
complex coefficients in the Laurent series 

I^ = ^ + ^ + ao, A^<4. (5.1) 

For the IR divergent triangles and boxes we use the analytic results of sec. For the UV 
divergent tadpoles and two-point functions we use the FF library of van Oldenborgh Q. 
For the finite integrals the coefficients a_2, a-i are equal to zero and we use the FF library 
for the coefficient uq. The code for the box integrals is unable to handle the phase space 
point where dety = and one is sitting exactly at the threshold given by the leading 
Landau singularity. The code assumes that all internal masses are real. An extension of 
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the code which also handles complex masses, appropriate for unstable particles, is a matter 
of analytic continuation and programming, rather than additional calculation. 

The code classifies the integrals in terms of the number of non-zero internal masses and 
then reduces the integrals to standard forms by relabelling the integral where necessary. 
Subsequently on the basis of the elements of the modified Cayley matrix the code identifies 
the appropriate divergent integral and evaluates it, or uses the FF library. 

To perform a numerical check for the divergent box integrals we make use the iden- 



tity |10| 



= ^ '^^3'' W + (3 - I))co/f , (5.2) 

where denotes the D dimensional triangle integral obtained from the box integral 
by removing the i-th propagator and the coefficients Cj are given by 

4 4 

a = ^{Y~^)ij , CO = ^ Ci . (5.3) 

i=i i=i 

The six-dimensional box, because it is finite, can be computed numerically using standard 
Feynman parameters. We introduce an prescription if needed and we set e equal to 
a small number. Assuming that the simpler, potentially divergent triangle integrals have 
been calculated correctly, using eq. (5.2) we obtain a rather powerful check of the numerical 



implementation of the box integrals, for both space-like and time-like values of the external 
invariants pf and Sij. 

The principal results presented in this paper, as well as the code can be downloaded 
from the website http : / / qcdloop . f nal . gov. 



6. Conclusions and outlook 

The essential new results of this paper are the classification of the infrared and collinear 
divergent triangle and box integrals, the calculation of the box integrals which were missing 
from the literature, and the provision of a code which returns a numerical answer for any 
one-loop scalar integral, divergent or finite, for four or less external legs. We believe that 
the problem of one-loop scalar integrals is now completely solved as far as next-to-leading 
order calculations are concerned. 

In conjunction with a procedure for determining the coefficients with which scalar 
integrals with four or less external legs appear in physical amplitudes we are in principle 
able to calculate the one-loop amplitude for any process. Amplitudes with massive internal 
lines, or massless internal lines, or both can now be treated in a seamless and uniform way. 
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A. Useful auxiliary integrals 

In this appendix we report on two integrals which were useful for the calculation of boxes 
11, 12 and 13. The first one is defined as 



Jo 



(m| -ml-pl + S23) 



li-mi - Pi) + (1 - l){ml - S23) 
' -7(1 - i)pI + iml + (1 - 7)ml 



X In 



) • (A.l) 



The result for the integral V is given in terms of the roots 7± of the quadratic equation 

-7(1 -7)pi + 7"ii + (1-7^1 = 0. (A.2) 

2 2 

Hence 7+7- = ^ and (1 - 7+)(l - 7_) = ^ and 

{pj + ml- ml) ± ^{pj + mi- mjf - Amjpi 
= ■ ^^-^^ 

Let us define the position of the pole in the integrand as 70, 

70 = o ' 0,0 • (A.4) 
^4 — •S23 — ^3 + P2 

In terms of these variables V is given by 

Viplphs23;mlml) = f ^7 ^Inf^p^W^) " ^^'^^ 

Jo 7-70 ^(l-7^)(l-7 

With this notation the result for V is 

T./ 2 2 2 2n t f-, 70 - 1 7^ \ T f. 70 - 1 7" \ 
V{P2,P3; S23; 'rn3,m4) = -L12 1 - L12 1 

\ 70 7^^ — 1/ V 70 7— 1/ 

+ (^) + Li. (^) + 2 Li. (i) . (A.6) 

If we set P2 = m| the pole is at 70 = 1 and we get 

V{mlpl-,S23;rnlml) = U2(^Y^^^ + Lb ( ) • (A.7) 

If we further set m| = ml = we have that 7+ + 7" = 1 and we obtain 

V(n^,A>.,;rrf.'^) = Li.(^) = -iln^ (^l^) . (A.8) 

In the limit — > we have 

-^^i + o{pI), -^-,!4 + <^(.pI)^ (A-9) 



7^4-1 734-1 
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and eq. (|A.6D simplifies to 

V{plQ- s,,-mlml) = Li2 (l - ^4^) + (l - ^ 



Li.(l- "f1-''^' ). (A.10) 

V m3(m|-S23)^ 



If we further set p\ = m| we get 



2 

y(mi 0; S23; mi, ml) = Li2 (l - ^) ■ (A.ll) 

Finally, if we also set m\ = mg = m? we obtain 

y(m^0;s23;m,^m2) = 0. (A.12) 

Now consider a related integral 

2 2 2 2n /"^ j (m| - ml - + S23) 
W^(P2'P3;s23;m3,m4) = / ^7— ^ — V— Ti TT^ V 

xlnf-^^^-^^^|+7"^3 + (l-7Kx 
V 7m| + (1 - 7)m| / 

Since 

W{pl,pl;s23]ml,mi) = V{pl,pl;s2'i\ml,ml) - y(p2, 0; S23; m|, m^) , (A.14) 

we find 

.rri 2 2 2 2n J. (. 70-1 \ / 70 - 1 ml 

b2,P3;S23;m3,m4) = L12 1 +L12 1 2 

V 70 / v 70 m^ 

-Li2fi- ^°-^ :\ )-Li2fi ^ 



70 7+ - 1 / V 70 7 

1 \ 1 \ m 



+ — j + Li4— j - Li^l - -Ij , (A.15) 
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